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Abstract 

The brane world scenario advocated by Arkani-Hamed et al. transmutes the hierarchy 
problem into explaining why extra dimensions have sizes much larger than the fundamen- 
tal scale. In this paper we discuss possible solutions to this problem by considering the 
compactified dimensions to be populated by a large number of branes in a crystal lattice. 
The experimental consequences of this scenario are described, including the presence of 
large energy gaps in the spectrum of Kaluza-Klein modes. 
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1. Introduction 

Around two years ago JjJ pointed out that it was consistent with known experiments 
for extra dimensions to exist with sizes of order a millimeter. The motivation for this 
observation came from string theory, where additional "curled up" dimensions are required 
for the consistency of the theory. String theory also allows for the existence of D-branes, 
which give a way to restrict the Standard Model fields to a three dimensional slice of 
the higher dimensional space. Without this additional entrapment of the Standard Model 
fields to a brane, large extra dimensions would be in immediate contradiction with atomic 
physics. 

The hierarchy problem becomes re-expressed as explaining why the size of the extra 
dimensions ro is much larger than the fundamental length scale. Denoting the fundamental 
scale by M* (which we can take to be of order 1 TeV), and the four-dimensional Planck 
scale Mpi = 10 19 GeV, one finds 

M 2 Pl = r%M:+ 2 , (1.1) 

for n flat extra dimensions with size tq. For n = 1, ro is required to be of cosmo logical 
scales, which is immediately ruled out. The n = 2 case requires ro of order one millimeter. 
There remains a hierarchy between ro and 1/M* ~ 10 _19 m. 

Current rounds of experiments , are beginning to place direct constraints on the 
simplest scenarios. Short-range gravity experiments now probe down to 200 //m strongly 
constraining the n = 2 case. Accelerator experiments constrain M* > 1 TeV ||. The 
strongest constraints arise from astrophysical considerations. Production of Kaluza-Klein 
modes in supernova SN 1987A places a bound M* > 30 TeV for n = 2 |[[. 

One of the main theoretical challenges in implementing some of the proposed large 
extra dimension scenario's is how to stabilize the size of the extra dimensions, without in- 
troducing additional fine tuning problems. This problem should be more readily addressed 
in the large extra dimension scenario, versus the traditional approach of compactification 
at the fundamental scale, because the analysis may be carried out at the classical level. 
A variety of solutions have been proposed in the warped compactification scenarios, 
where by warped we mean not only the Randall-Sundrum warped metric scenarios []T0| |13 1 



but also cases where some other field, eg., a scalar as in [§, has a non-trivial profile in the 
extra dimensions. 
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Another proposed solution in the unwarped compactification context of n is to have a 
large number of branes iV (one of which would be our world) interacting in such a way that 
they form a crystal lattice in the internal dimensions |14j . The interbrane separation could 
be the size of the fundamental length scale, but the size of the extra dimensions would 
nevertheless be large enough to solve the hierarchy problem for a large enough number of 
branes. While it is still not clear how to realize such a brane lattice crystal from a more 
fundamental theory such as string theory, these models are nevertheless of interest from 
a phenomenological point of view, especially in light of the fact that they could be tested 
experimentally in the near future. 

One is still left with the problem of explaining the large integer N. We take the point 
of view that replacing a fine tuned continuous parameter by a large integer parameter is 
an improvement, as one can set the integer iV once and for all using initial conditions, and 
under suitable circumstances, this integer will be stable with respect to time evolution. 
Since the hierarchy is set by a conserved number N, it is automatically stable with respect 
to radiative corrections. 

In the following we begin by discussing the approach to brane crystals of Arkani- 
Hamed et al. ]14| . They propose a number of scenarios for stabilizing the extra dimensions. 
One particularly compelling example does this without introducing additional fine tuning, 
aside from the usual problem with the four-dimensional cosmological constant. A problem 
with this scenario is the presence of unbalanced charge on a compact space. Balancing the 
charge leads us to consider a scenario where neutral non-BPS branes interact via a nearest 
neighbor potential. We show this does lead to a natural solution of the hierarchy problem. 
The crystal potential leads to a distinctive experimental signature for this scenario - the 
existence of large energy gaps in the Kaluza-Klein spectrum. Such energy gaps have also 



been noticed in the Randall- Sundrum scenario brane lattices discussed in [15,16 



2. Brane Crystal Review 



We begin by reviewing the model considered in |14j| . They consider a 3-brane em- 
bedded in a universe with 3 large spatial dimensions and n small spatial dimensions. The 
system is described by a bulk action 



'bulk 



J d 4+n x^/-detg 4+n {Ml +n TZ + A - C matter + • • •) (2.1) 



and a brane action 



Sbrane =~J <?X ^~ det g 4 (/ 4 + •••), (2.2) 



where Cmatter is the Lagrangian of the bulk matter fields and the ellipses denote higher 
derivative terms which may be dropped at low enough energies. Here g<± denotes the 
induced metric on the brane, g^+n denotes the bulk metric, A is the bulk cosmological 
constant and f 4 is the brane tension. Interaction terms between the branes are not in- 
cluded. In the next section we discuss scenarios where such interaction terms are relevant. 
The metric is assumed to take the form 



ds' 



(— (dt 2 - R 2 g l3 dx l dx j ) - r l g IJ dx I dx J , (2.3) 



where R = R(t) is the scale factor of the three large dimensions and r = r(t) that of the n 
small dimensions, with tq = r (0). 0. Also we use lower case Latin indices • • • to denote 
the 3 large spatial dimensions and upper case Latin indices I, J, ■ ■ ■ for the n small spatial 
dimensions. 

Inserting this form of the metric into the bulk ( [2.1| ), and brane ( |2.2j ) actions and 
integrating over the spatial coordinates results in 

S = I dtR> | M-V; I -6 m + 2S=±2 (t) 2 )- C-l) 2 " Vtot (r) I , (2.4) 




after an integration by parts. In obtaining this result we have assumed that J d 3 x^det gu = 
1 and J d n x^Jdetg~ij = 1. The potential V to t(r) is given by 

V tot {r) = Ar n - Kn(n - l)M, n+2 r n " 2 + f , (2.5) 

where the term proportional to k arises from the curvature of the n dimensional space. 
For an n-sphere k = 1 and for an n-torus k = 0. A similar term could be added for the 
three large spatial dimensions, but for large R would be negligible so we drop it. 
For static solutions the equations of motion are simply given by 

Vtotiro) = 0, VUro) = 0, (2.6) 



1 [14] do not include the conformal factor r~ n multiplying the large dimensions. We find it 
more convenient however to include it as it will lead to diagonalized kinetic terms for r and R in 
the action below and is the ordinary conformal factor appearing in Kaluza-Klein reductions. 
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assuming that R(t) = Ro is constant. Note that with the potential (|2.5|) r is fixed entirely 
in terms of the bulk interaction. For the above potential (|2.5| ) these equations constrain A 
assuming that ro is chosen to solve the hierarchy problem (|1 . 1| ) . Specifically let's assume 
that the brane tension is set by the higher dimensional fundamental scale M* so that 



JU-n+2 / M \ 4/n 



N 



f 4 « M 4 . From (|2.6|) , assuming the internal space is an n-sphere, it follows that 

^1 „ M n + 4 (M*_ 

rg * \M Pl 

where we have used (|T71|). 

For A" branes with equal tensions, the / 4 term in (|2.5| ) is replaced by Nf A . Solving 
the V to t(r ) — equation leads to 

M \ 2(n-2)/n 

. (2.8) 
M*J K ! 

This varies from 1 for n = 2 to 10 20 for n = 6. Note one is still left with an extra fine 
tuning problem, in order that the bulk cosmological constant satisfy the relation (|2.7|) . Also 
note the brane number A" plays no role in fixing the size of the extra dimensions - this is 
entirely determined by fine tuning A. A" is fixed only by requiring the four dimensional 
cosmological constant vanish. 

For a single brane, stability of this solution follows from the condition V^ t (ro) > 
where ro solves the equations of motion (^]6|) . This is straightforward to see by expanding 
the action to quadratic order in the small perturbation 5r where r = ro + 5r and then 
rescaling 5r so that the kinetic term is canonically normalized. The mass can then be read 
from the action and is given by 



< = -7- I <^ " 72 J ■ (2-9) 



1 rlV[' ot {r Q ) 
n(n + 2) M pi 

Evaluating this for ( |2.5| ) gives m r w 1/tq. Note the static equations of motion 
and stability condition do not depend on the specific form of Vtotix), and yield strong 
constraints on the parameters of more general potentials. 

This analysis presumes that derivative couplings of the radion to higher spin Kaluza- 
Klein modes may be neglected, which is not true in general. However this should not 
change the qualitative conclusions. The analysis also neglects the Hamiltonian constraint, 
which would take the form 

m." + V;. | -e ( £ ) 2 + (L) 2 ) + Vtot(r) = o . (2 , ) 




This tells us if we really considered perturbations independent of the three noncompact 
spatial dimensions, we would generate an non-zero energy density everywhere in space, 
which would lead to expansion or contraction of R. This is easily remedied by generalizing 
to perturbations localized in the noncompact spatial dimensions. 



3. Interbrane Forces 

We now consider in more detail the effect of an interbrane potential on the above 
analysis. We continue to work in an approximation where the compactification is not 
warped, i.e. the 4 + n-dimensional metric does not depend on the internal coordinates 
x 1 . This presumes the brane separation will be stabilized at a size parametrically larger 
than the fundamental length scale 1/M*, so that treating gravity at the classical level is 
sufficient. We also continue to treat perturbations homogeneous in the spatial directions, 
with only time dependence, with the metric ansatz ( |2.3| ). 

Let us generalize the brane action to N branes, with a Born-Infeld type action 



(fc) 



Sbrane = " / ^xJdet - 9l jd ^X^xQ (/ fc 4 + ^ V brane (r\X® ~ X 

k=l \ Kk 

(3.1) 

with an interbrane potential Vb rane . Here we have assumed the metric is g = g(t) induced 
on each brane is the same, and we use X^ = Xr^ (t) to denote worldvolume fields corre- 
sponding to the brane positions in the transverse space. Indices fj,, v label the coordinates 
(t, x l ). The brane potential has been chosen to depend on the proper distance separating 
the branes. This leads to the factor of r in potential term. In writing the Born-Infeld action 
we have used the worldvolume diffeomorphism invariance to fix the gauge X? k ~. = x M . Note 
that had we included spatial dependence of the fields, the brane potential would take a 
complicated non-local form, that is difficult to write down explicitly. A useful analogy for 
the branes interacting with the bulk gravitational field is gravitational waves interacting 
with a resonant gravitational wave detector, |17[] , as we will comment further below. 

A natural candidate for Vj, rane is a simple Coulomb coupling. This leads to the most 
interesting brane crystal scenario studied in [nj] with "non-extensive bulk cosmological 
constant", where the hierarchy problem was solved without the additional fine tuning 
associated with the bulk cosmological constant. The radion was stabilized in the infra-red 
using a negative curvature term in the internal space, and the Coulomb force was used 
to provide a short distance stabilizing force. The difficulty with this picture is that a 
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collection of like charged branes on a compact space carries infinite vacuum energy, since 
the electric flux has nowhere to end. 

To remedy this, one could consider brane configurations with zero charge per unit 
cell. Of course, once branes with opposite sign are present there will be attractive forces. 
For supersymmetric D-branes in string theory, oppositely charged branes will annihilate. 
Furthermore, to our knowledge, there are no known neutral and stable branes. 



The description of D-brane charges as living in K-theory groups [T^] however predicts 
the existence of stable non-BPS branes which carry charge in a finite, or torsion, group. 
This charge is not associated with a gauge symmetry and there is therefore no Gauss law 
preventing us from considering iV such branes on a compact space. It remains an open 
question as to whether such branes could be used to construct a stable lattice configuration. 

For the moment we take a phenomenological point of view and assume that a stable 
lattice can be constructed. These objects will then not experience a Coulomb interaction. 
A Van der Waals interaction is one natural interaction between such objects, arising from 
the interaction of induced electric dipole moments, falling off like l/r 2n . To obtain a 
stabilizing potential, this must be combined with a hard-core repulsive interaction. Taking 
our motivation from molecular crystals, a possible potential would be the n-dimensional 
version of the Lennard- Jones potential 

x) ~ (x) ' (3-2) 

where f3 and 7 are both 0{1). One could also imagine an ionic lattice of branes, with 
screened Coulomb interactions. An importance difference with the Coulomb force example 
is that now V^rane will scale like the number of branes N, rather than iV 2 since nearest 
neighbor interactions will be dominant. The precise form of the potential will not be 
important for what follows. 

We now want to show that the size of the internal space is fixed in terms of the number 
of branes N, rather than by using the bulk quantities k and A, which generally introduce 
extra fine tuning problems. We shall therefore set k and A to zero, the former implying 
that the compact extra dimensions are flat and for simplicity we take their geometry to be 
the torus (Si) n . 

Consider an interbrane potential of the form 

V hrane {r\X {k) -X®\) = M$v{\XW -X®\rM.) , (3.3) 
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where v(x) is not fine tuned. The fundamental scale M* sets the scale of the interaction. 
We assume v(x) is short ranged, so only nearest neighbor interactions are dominant. The 
static equations of motion are then roughly given by v(xq) = 0(1) and v'(xq) = 0. These 
have solutions xq = 0(1), or more explicitly 

Ax « ~tf- > (3-4) 

where AX denotes the coordinate separation between nearest neighbor branes and a is 
a constant of order 10 or so, which we will discuss momentarily. Using the fact that the 
coordinate periodicity around any of the S 1 ^ of our extra dimensions is 1 and summing 
up the AJ's along one of these dimensions yields the static value of r, 

N 1 / n 

ro - . (3.5) 

We find therefore that the size of the internal dimensions is set by the number of branes 
along with the fundamental scale. The value of ro however was already fixed in (|T7l|) and 
thus yields the necessary number of branes 

N „±(«p)' = (3.6) 

a = 1 corresponds to one brane per fundamental volume in the internal space, saturating 
the number of branes. We have therefore required that a is of order 10 or larger (but not 
so large that we have another fine tuning problem) in order that classical gravity be a good 
approximation. 

The effective potential induced for r takes the form 

V tot (r) = Nf 4 + NVtraneir/N 1 ^) , (3.7) 

where integrating out the brane coordinates sets the coordinate distance between neigh- 
boring branes to l/N l / n . Vanishing of the four-dimensional cosmological constant requires 
one fine tuning, corresponding to V to t = 0, but note no additional fine tuning is needed. 
The mass of the radion may be obtained using (|2.9|), which gives 



m r = Mja n/2 . (3.8) 

We have so far taken vanishing bulk cosmological constant. This is expected at tree- 
level if supersymmetry is unbroken in the bulk. However if supersymmetry is broken on 
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the branes then a cosmological constant will appear at one loop and could change the 
results. If the breaking takes place at the fundamental scale M* on the branes then the 
induced mass splittings in the bulk are given by a tree-level gravitational effect 



Am 2 » N-S- = — , (3.9) 

where the last expression was obtained by evaluating tq at ( |3.5| ). By dimensional analysis, 
this induces a cosmological constant ^quantum ~ (Am 2 )( 4+n )/ 2 so that the potential gets 
a contribution of 



M 4+n / Tn \ n 

V quantum (r) - MtN {-jz^rjA , (3.10) 



where xq was defined above. Therefore at xq ~ a the induced cosmological constant 
contribution to the potential is subdominant and our original estimates above still apply. 
We expect the brane crystal scenario will only work when the number of extra di- 



mensions n > 3. This follows from [19] where it is shown that classical crystal lattices 
do not exhibit long-range order in dimensions two or less. For M* w 1 TeV, (|1 . 1|) places 
vq « 10 -7 m for n = 3. 

4. Experimental Consequences 

To understand the experimental consequences of our scenario we must investigate the 
spectrum of the theory as well as the couplings to Standard Model fields. To do this it is 
convenient to first fix the (4 + n) dimensional diffeomorphism invariance. For linearized 
perturbations Hmn about a flat metric ()mn (with M,N labeling the 4 + n-dimensional 
space) an infinitesimal diffeomorphism generated by a vector £m acts on the metric as 

hMN — > hMN + VmCjv + Vat^m , (4.1) 

where we have decomposed the metric into a background part ()mn and a fluctuation hMN 

as 

9mn = 9mn + hMN ■ (4.2) 



On the -X"(fc)'s this will act as 



X^X 1 -^. (4.3) 



Part of the gauge symmetry can be fixed by demanding 



d M h MN = , (4.4) 

where 

h>MN = h-MN — ~9MNh , (4.5) 

and indices are being raised and lowered with the background metric. This gauge choice 
does not fix the diffeomorphism invariance corresponding to vectors satisfying the (4 + n)- 
dimensional wave equation d 8m£,n = 0. 

We begin by expanding the Born-Infeld piece of the brane action for a general pertur- 
bation 

J (fxy/detG^ = J <i 4 x v /det^(l + ^g^n^ 
where is the full induced metric on the brane, and 



(4.6) 



, dX 1 , dX 1 _ dX 1 dX J , if . 

Since picks up a nontrivial second-order contribution under diffeomorphisms of the 
form ^ = 0, ^ 

V V + Ki 9 ^ 1 + huid^ 1 , (4.8) 

as may be seen by expanding (^4.1|) , it is convenient to write h^ v = j M „ + h^ih^, so will 
then be invariant under such diffeomorphisms. This redefinition also makes clear that (|4.6| ) 
gives rise to a mass term for the four-dimensional vector fields h^j . A different approach to 
seeing the vector fields become massive is discussed in |1|] . However since we are tuning the 
four-dimensional cosmological constant to zero, the overall coefficient of Q4.6| ) will vanish, 
and the vectors will be massless (as is the graviton). 

To examine terms arising from the brane potential, we again restrict to perturbations 
independent of the spatial directions. We gauge fix the fluctuations of the brane coordi- 
nates Xi-, to zero, as explained in more detail in appendix A. At quadratic order in the 
fluctuations, then only the radion modes hjj couple to the brane coordinates XL, in the 
action. 

The modes independent of the internal dimensions will therefore be a massless graviton 
j^v, a set of massless vector fields hi, and a set of massive radion scalars hjj with mass 



given by ( j3.8| ), following through the same calculation. Each of these modes will be at 
the bottom of a tower of Kaluza-Klein states which are standing waves in the internal 
dimensions. These may be treated in the same way as Bloch waves |f2Q| . For the low 



lying modes, the effects of interactions may be neglected (at least for sufficiently large 
a). This is precisely analogous to the case of gravitational waves propagating through a 
resonant detector, where one needs to go to next order in the equations of motion to see 
the effect of gravitational interactions on the response of the detector [|I7] . This gives rise 
to a typical Kaluza-Klein spectrum for the spin- two and vector modes rrik = |/c/ r o| ? while 
for the radion modes m 2 = m 2 + (k/ro) 2 . 

It is also interesting to calculate the energy band gap at the edge of the Brillouin zone 
boundary, where we have standing waves commensurate with the lattice spacing of the 
crystal. In general, interaction effects will become large there. For the radion modes, we 
can get a reasonable estimate of this band gap by taking into account only the interaction 
of the radion through the brane potential. For plane wave modes propagating in the J'th 
direction, the equations of motion are the same as that of an electron moving in a periodic 
array of delta function potentials. This is a special limit of the Kronig-Penney model |2lj] . 
The wavef unction takes the form 



r(x) = e lkx u(x) , (4.9) 

where x is the x 1 in question, and u(x) is periodic under lattice translations. Solving the 
equations of motion for u(x) a linear combination of e ±lKx yields 



1 sm(Ka) . ir 

cos(ka) = cos(Ka) + ^ , (4.10) 



where a is the lattice spacing tq/N 1 / 71 . Here Ka is to be identified with ctm/M*, where 
m is the mass of the mode. For the first Brillouin zone boundary k = re/a. Solving this 
equation yields the band gap Am 2 ~ M 2 /a n . This is the same form as the mass gap of 
the radion near k = 0. 

Higher order interactions also lead to a band gap for the spin- two and vector modes 
at the Brillouin zone boundary. Estimating the energy difference between a standing wave 
with nodes on the branes versus a standing wave with peaks on the branes leads to the 
same calculation as in (jO|). We therefore expect the band gap Am 2 to be of roughly the 

10 




ka 

Fig. 1: Illustration of expected band structure in the spectrum of Kaluza- 
Klein states. The mass squared is plotted versus the wavevector in the 
internal dimension k times the lattice spacing a ~ a/M*. The orders of 
magnitude of the scales of the band width and band gaps are indicated. 



same order of magnitude as for the radion modes. The picture of the band structure that 
emerges is illustrated in figure 1. 

The Standard Model fields are coupled to the induced metric on the brane via the 
usual covariant couplings. Expanding these terms about the background metric we find 
that the coupling between the bulk metric fluctuations discussed above and the Standard 
model fields will be suppressed by 1 /Mpi. The analysis of the phenomenological constraints 
of [0 will therefore carry over to the brane crystal model unchanged. 

Acknowledgments We thank A. Houghton, R. Myers, R. Pelcovits and L. Ran- 
dall for helpful discussions. This research is supported in part by DOE grant DE- 
FE0291ER40688-Task A. 
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Appendix A. Gauge fixing the brane coordinates 

As noted above, the linearized Einstein equations in covariant gauge d M JiMN have 
a residual gauge freedom associated with diffeomorphisms satisfying d Om£,n = 0. This 
freedom can be used to set the fluctuations of the transverse brane coordinates to zero given 
some assumptions. Specifically the brane coordinates transform under a diffeomorphism 

as 

Xf fc) (s") ^X( k) (x^-e(x^x j ) , (Al) 

where we have assumed that £ M = 0. To set the brane coordinate fluctuation to zero 
therefore only requires fixing £ J on the brane, and it's extension off the brane can be chosen 
at our convenience. In this case we choose it's extension off the brane so that it solves 
the (4 + n) -dimensional wave equation. In general this will not be possible. This follows 
by noting that the £ 7 live in the bulk spacetime and therefore must satisfy the periodicity 
conditions of the extra dimensions. As a result the £ 1 cannot have arbitrary dependence on 
the x M coordinates if they are to satisfy the wave equation. To be explicit consider a mode 
decomposition of £ J where one of the modes satisfies d^d^ 1 = m 2 ^ 1 . For £ 7 to satisfy the 
(4 + n)-dimensional wave equation we then need d J dj^ 1 = m 2 ^ 1 . For generic values of m 
there will be no solution to this equation that satisfies the required periodicity conditions 
in the extra dimensions. Rather the 4-dimensional mass m will be quantized according 
to m = \J q 2 + • • • + Qn/ r o f° r arbitrary integers q^. There is no such constraint on the 
dependence of the brane coordinates on however because they are functions only of the 
x M . We assume for simplicity nevertheless that the brane coordinate fluctuations can be 
gauged away in this manner. In other words we assume that the mode decomposition of 
the brane coordinates contains only fluctuations with the quantized masses given above 
therefore allowing us to gauge them away. 

Appendix B. Alternate derivation of linearized equations 

Given these gauge conditions described in appendix A, we now give a more general 
derivation of the linearized equations of motion describing the metric fluctuations. We 
could try to write down a potential term describing the interactions between the branes 
and then gauge fixing as described above, but this turns out to be somewhat subtle. We 
therefore argue using symmetry considerations. The bulk contribution comes only from 
the Einstein-Hilbert term in the action (|2.2| ). In the gauge ( |4.4| ) it is well known that this 
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contributes only d p dpliMN to the linearized equations of motion. From the brane terms 
in the action we expect the linearized equations of motion to contain a sum of 5-function 
terms in the extra dimensions corresponding to the fixed brane positions with coefficients 
determined by symmetry and dimensional analysis. Specifically we find 

M 4 

d M d M h^ = T7T V + b^-g^h^iX 1 - X( k) ) 

M 4 

d M d M h,j = -f ^c W V^(^ -X( k) ) (B.l) 
M 4 

d M d M h u = jjt Y,( d (k)hiJ + e^guhn^Hx 1 - X[ k) ) 
Pl (fc) 

where = g^ u h^ u , h n = g IJ hu, and the coefficients ar^, e^) are all 0(1). The 
overall factor of M 4 /Mp t on the right-hand-sides of all equations is easy to understand by 
going to a coordinate system X 1 = r^X 1 . In these coordinates M* is the only dimensionful 
parameter so that M^~ n would have to be the overall coefficient following from dimensional 
analysis. Going back to the X 1 then yields the above coefficient. 

Decomposing the fluctuations Jimn into eigenstates of the 4-dimensional wave operator 
yields 

d M d M h MN = (-ro 2 did T + m 2 )h M N (B.2) 

for a mode with 4-dimensional mass m. It is now straightforward to estimate the spectrum 
of metric fluctuations following the discussion around ( |4.9| ). In particular taking a plane 
wave ansatz for the metric fluctuations one recovers the relation ( |4.10| ) from which we find 
that the first excited state and mass gap energies will be or order m 2 ,Am 2 w M%/a n 
respectively. 

The discussion so far applies for all three equations of motion in (|B.1| ), so in particular 
it implies that the lowest energy fluctuation of would have four-dimensional mass of 
order M*/cW 2 . This is of course unacceptable if we are to recover Newtonian gravity on 
our brane. For the equation therefore we must tune the ar^) and bn,} coefficients so that 
we have a massless fluctuation, or massless four-dimensional graviton. This corresponds 
to fine tuning the four-dimensional cosmo logical constant to zero. Note however that this 
fine tuning will not in general imply that the low lying vector fluctuations h^j or radion 
fluctuations hu will be massless. 
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